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Two, the most simple cases of special-relativistic flows of a viscous, incompressible fluid are 
considered: plane Couette flow and plane Poiseuille flow. Considering only the regular motion of 
the fluid we found the distribution of velocity in the fluid (velocity profiles) and the friction force, 
acting on immovable wall. The results are expressed through simple analytical functions for the 
Couette flow, while for the Poiseiulle flow they are expressed by higher transcendental functions 
(Jacobi's elliptic functions). 
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. A recent development of theoretical astrophysics revealed a wide class of hydrodynamic and hydromagnetic flows 
7-H ' having the close relevance to a number of astrophysical objects such as stellar winds H], accretion flows 0], outflows 
of matter in Active Galactic Nuclei (jets) Q etc. The spatial symmetry of these flows is different, thermodynamic 
state is sometimes unusual and velocities are often relativistic. Relativity definitely influences the physical processes 
that take place in such flows and leads to various outstanding appearances of these objects. 

The subject of the present paper is to study how the relativistic velocities affect physical conditions and dynamics 
of the most simple model hydrodynamic flows. For the present consideration we chose such examples of flows, which 
■ being simple enough are, at the same time, relevant to the most astrophysical cases due to their simple geometry. 

In particular, we are going to study the special-relativistic generalization of a couple of examples of the well-known 
classic hydrodynamic flows M: plane Couette flow and plane Poiseuille flow. We consider only regular motion in these 
flows. In both cases we find exact solutions of the problem. 

In particular, for the plane Couette flow we determined the velocity distribution throughout the flow. The velocity 
i 1 profile is expressed analytically and appears to be nonlinear. The nonlinearity noticeably increases with the increase 
JL| | of the velocity of the moving plane and becomes almost "steplike" for vo— 1. We have found the average velocity 
of the flow that in nonrelativistic limit tends to vq/2 (as it should be), while for vq~1 it tends to unity. We have 
also determined the friction force acting on the immovable wall, which is exactly equal to the corresponding classic 
i— I . (nonrelativistic) value. 

For the plane Poiseiulle flow the velocity distribution was derived by the expressing of the solution of the correspond- 
^ , ing differential equation through the higher transcendental functions. In particular, the velocity profile is expressed 
by Jacobi's elliptical functions. The general property of this solution is that with the increase of the pressure gradient 
value (the very gradient that induces the whole flow) the profile becomes more and more "flattened." 

In the concluding section we discuss the obtained results in the context of their possible relevance to the future 
modeling of real relativistic astrophysical flows. 
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II. MAIN CONSIDERATION 



Henceforth, we shall use the following notations: greek indices will represent spacetime components, while latin 
ones will be used for spatial components. We shall use geometrical units, so that G = c = 1. The signature of the 
spacetime (Minkowskian) metric: 

ds 2 = -dr 2 = -dt 2 + dx 2 + dy 2 + dz 2 , (2.1) 

is chosen as (— + + +)■ 

We shall consider motion of the incompressible, viscous fluid with the following stress-energy tensor ||: 

T a/3 _ gJjajj/3 + p^fi _ ^^(i _ (2.2) 
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where all notations are standard |5[. In particular, e is the density of the mass-energy, P — the pressure, U a is the 
4- velocity field of the flow, h a @ is the projection tensor, defined as: 

h c*0= v *0 + u a uP, (2.3a) 

r] and £ are the coefficients of the shear and bulk viscosity respectively. 9 is the "expansion" of the fluid world lines: 

9=U£, (2.36) 

and a a P is the shear tensor: 

- \eh a p. (2.3c) 

For the fluid with nonrelativistic temperature, density of the mass-energy e is equal to the "usual" density of 
matter p in the proper frame of reference. Since a fluid is thought to be incompressible: p = const and the continuity 
equation: (pU a ). a = 0, clearly leads to 9 = ("expansion" of the incompressible fluid is equal to zero). 

Under mentioned circumstances the stress-energy conservation law may be generally written as: 

T a % = 0. (2.4) 



A. Plane Couette flow 



Let us consider the plane flow of the incompressible viscous fluid, which is situated between two infinite parallel 
planes one of which is moving with arbitrary velocity vq. Let an X axis be parallel to the planes and an Y axis 
be normal to them so that y — in the plane at rest, while y = L for the plane in motion. The symmetry of the 
problem implies that the sole nonzero component of the fluid velocity vector field v{f, t) is v x — const(x). Naturally, 
the 4-velocity nonzero components should be equal to: 

ff* = (1 - t£)-V2= 7 , (2.5a) 



U x = jv x . (2.56) 

Taking into account (2.5) one easily finds that the x-th, nonzero component of the stress-energy conservation equation 
(2.4) reduces to: 

and its general solution, evidently, is: U x = Ay + B where A and B are some constants of integration. Since the fluid 
at the immovable plane (y — 0) should be at rest, B = and then from (2.5b) we can obtain expression for the x-th 
component of the 3-velocity, which is: 



Ay 



(2.7) 



Another boundary condition: v x (y = L) = vq helps to find the remaining integration constant A: 

A = Jovo/L, (2.8) 

and (2.7) may be written explicitly as: 
as for 7(2/) we get: 

7(y) = v / i+PF7- (2 ' 10) 
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Now, let us calculate an average velocity of this flow v. According to the definition we should have: 

v=j J v x (y)dy 

Taking into account (2.9) and performing necessary integration we get the following result: 



VO Vq 

When vq<^1 (2.11) leads to the self-evident asymptotic result v~vq/2 [Q. For the case vq— >1 the average velocity 
V— >1 as it, certainly, should be. 

It must be noted that the tangential friction force, acting on the immovable (y = 0) plane, defined as: 

f=2n\a xy \ = ri\dU x /dy\ = nv /L, (2.12) 

coincides with the classic, nonrelativistic result Q. 

B. Plane Poiseuille flow 

Let us consider, now, a flow of the viscous, incompressible fluid between two parallel infinite immovable planes, 
induced by the existence of the pressure gradient along the flow axis. Let X axis be coincided with the symmetry 
axis of the flow. Therefore, the planes are situated at y — ±L. The symmetry of the problem implies that: 

U a = [U t ,U x =U = f(y),Uy = 0]. (2.13) 

Note that the y-th component of the stress-energy conservation equation leads to P = const{y). The x-th component 
of the same equation reduces to the following expression: 



DP ^d 2 U 
dx dy 2 



^ + U 2 )— -to =0- (2.14) 



As far as U = const{x) the gradient of the pressure along the flow should be constant: dP/dx — const. Furthermore, 
this constant should be negative if U is taken to be positive. Expressing the second derivative appearing in (2.14) as: 



d 2 U dU d fdU 



dy 2 dy dU \ dy 



we can rewrite the equation as: 
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V 


dx 



(2.15) 



{l + U 2 ). (2.16) 



Introducing new dimensionless variable: e=(2\dP/dx\/3r]) 1 ^ 2 y we reduce (2.16) to the following first order differ- 
ential equation: 

dU 



r ^-(U 3 + 3U + const), (2.17a) 

where const is some constant of integration. It may be determined through the obvious condition, arising from the flow 
symmetry: (dU/de) £ =o = and U{e = 0)=U ma x- These conditions together imply that const = —U m ax 3 — 3[/ maa: . 
Thus (2.17a) may also be written as: 

dU 



;)e = y/(U max - UWLx + UmaxU + U 2 + 3), (2.176) 

This equation is solved in higher transcendental functions — Jacobi's elliptical functions (see, for details || and 
Appendix of this paper). Taking into account that according to (A. 5): 

-Ae = F(tp\m), 
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where 



A 2 = v/3(l + O 



a/3 

1 + — y m aa; 



(2.20a) 
(2.206) 



COS<y9 = 



X 2 — {Umax ~ f ) 



A 2 + (f/„ 



for U(e) we find the following expression: 



U(e) = U rnax - 



Xsn(Xe) 



n 2 



(2.20c) 



(2.21) 



1 + cn(Xe) 

As for the x-th component of 3- velocity V=v x (y) we can simply derive it through (2.20) via the obvious relation 

U(s) 



V(s) 



(2.22) 



We see that with increasing of V max velocity profile becomes more and more " flattened out" . The similar tendency 
is clearly seen also for Couette flow, considered in the previous section. 



III. CONCLUSION 



In this study we consider two kinds of the special-relativistic flows of incompressible, viscous fluid. In particular, 
we examine the most simple two-dimensional (plane) model flows between two parallel infinite planes: Couette flow 
and Poiseuille flow. In both cases we find distribution of the velocity throughout the flow (velocity profiles) and in 
the former case we have also found the average velocity of the flow and the friction force, acting on the immovable 
wall. The obtained solutions are exact: for Couette flow it is analytic, while for the Poiseuille flow it is expressed by 
Jacobi's elliptic functions. 

As the further step of the investigation one may consider a study of hydrodynamic instability of these flows. It is 
presumable that the relativity may introduce new qualitative features in the instability dynamics as well. However, 
such a study is beyond the scope of the present paper and will be considered elsewhere. 

Since the geometrical properties of considered model flows are the most simple we think that the results obtained 
in the present study and the mathematical methods of research that we are here developing, may be relevant and 
useful in some kinds of prototype astrophysical relativistic flows with comparable properties. Here we mean the rela- 
tivistic regions of accretion discs of compact objects, accretion flows and stellar winds with spherical or quasispherical 
symmetry, outflows of matter in Active Galactic Nuclei, etc. Certainly, in some of these real astrophysical flows 
special-relativistic effects will be intricately interlaced with general-relativistic effects evoked by strong gravitational 
fields and electromagnetic effects (for plasma flows) arising due to the existence of superstrong magnetic fields. We 
hope that the present study may be of some use as the standing point in the future investigations of these complicated 
types of relativistic astrophysical flows. 



IV. ACKNOWLEDGEMENTS 

My research was supported, in part, by International Science Foundation (ISF) long-term research grant RVO 300. 

V. APPENDIX 

In the general theory of Jacobi's elliptical functions it is known that the integral: 
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(where P(£)=i 3 + ait 2 + a^t + 03 is the third order polynomial of t, having only one real root t = (3) may be reduced 
to the Jacobi's elliptical integral of the first kind. For this purpose one must introduce the following coefficients: 



X 2 =y/F^), (A3a) 



m=sin 2 a = \- \^ML, (AM) 
2 8v^W 



and use, also, the following substitution: 

A 2 + (fi - t) ■ 



>?-(P-t) .... 
coscp = y — -(. (A A) 



In these notations the initial integral may be written as: 

dt 1 r d6 



(A5) 



The latter integral is equal to the Jacobi's elliptical integral of the first kind: F((p\m). 

In our case a\ — 3; 02 = and 03 = const. Despite that 03 is an arbitrary parameter, the sole real root of the 
polynomial appearing in (2.17) should be the value of U(e) on the flow axis, since the symmetry of the flow implies 
that at e = 0, dU/de = and at this point U — U„ 



' max- 
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